Prospects for observations of high-energy cosmic tau neutrinos 
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O ■ We study prospects for the observations of high-energy cosmic tau 

neutrinos (E > 10 6 GeV) originating from proton acceleration in the 
cores of active galactic nuclei. We consider the possibility that vacuum 
flavor neutrino oscillations induce a tau to muon neutrino flux ratio 
CC) | greatly exceeding the rather small value expected from intrinsic produc- 

tion. The criterias and event rates for under water/ice light Cerenkov 
neutrino telescopes are given by considering the possible detection of 
downgoing high-energy cosmic tau neutrinos through characteristic dou- 
! ble shower events. 

^ ■ PACS number(s): 95.55.Vj, 14.60.Pq, 13.15.+g, 98.54.Cm 
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I. INTRODUCTION 



OS. 

Neutrino astronomy is now an emerging field and entails the need to have improved flux 



estimates as well as a good understanding of relevant detector capabilities for all flavor neu- 
trinos, particularly in light of the recent growing experimental support for flavor oscillations 
Several high-energy cosmic neutrino (>~ 10 6 GeV) detectors based on under water /ice 
muon detection are now at proposal or construction stages [fj] and alternative techniques 
for 10 9 GeV neutrinos are being considered through coherent radio or acoustic M pulses 
as well as through horizontal air shower measurements with conventional arrays || or with 
fluorescent light either from the ground |J[7| or from orbiting detectors ||. A number of as- 
trophysical high-energy neutrino sources, such as Active Galactic Nuclei (AGN), have been 
discussed in the literature and predicted to produce fluxes that could be detected in some 
of these detectors ||. In the event of a successful detection of high-energy astrophysical 
neutrinos the range of parameter space for flavor oscillations that can be tested could be 
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considerably enhanced provided a flavor identification can be done. 

Cosmic tau neutrinos are possibly the easiest flavor to identify above 10 6 GeV. Two 
ideas have already been put forward based on the short decay lifetime of the r produced in 
charged current interactions. There is a suggestion of measuring 10 6 GeV v T flux through 
double shower (double bang) events in under water/ice Cerenkov telescopes flT0| . A more 
recent suggestion is to detect a small pile up of upgoing /i-like events in the 10 4 — 10 5 GeV 
range with a fairly flat zenith angle dependence [IT|. On the other hand the LPM effect, that 
lengthens electromagnetic showers in water and ice, has been suggested to separate electron 
neutrino charged current interactions from the rest. This could be done with Cerenkov light 



detectors or with the radio technique for energies above 2 • 10 7 GeV ||12||. It is conceivable 



that a combination of several of such techniques will allow the establishment of neutrino 
flavor ratios at energies above 10 6 GeV. 

We will concentrate on cosmic tau neutrino detection in this article. Specifically, we 
discuss in some detail, the prospects for detection of high-energy cosmic tau neutrinos orig- 
inating from the cores of AGN. We consider vacuum flavor oscillations as an example to 
illustrate the possibility offered by detectors in construction to distinguish between differ- 
ent neutrino flavors. For absolute event rates we use upper limit flux calculations as an 
example |13| and consider 1 km 2 detector sizes which are now being planned 0. Both the 
energy ranges of interest and the relative numbers of r- and /x-like event rates are however 
independent of the assumed normalization of the neutrino fluxes. 

We show that for the chosen neutrino flux, a km 2 size surface area under ice/water 
Cerenkov light neutrino detector may be able to either set useful upper limits or may obtain 
first examples of the high-energy tau neutrinos originating from this cosmologically distant 
astrophysical source. The plan of the paper is as follows: In Section II, after a brief discussion 
of intrinsic production mechanisms of high-energy muon and tau neutrinos in AGNs, we 
estimate the relevant vacuum flavor oscillation probability. In Section III, we discuss in 
some detail the detection technique making use of the double shower structure of the tau 
neutrino charged current interactions and calculate the expected event rates for a typical 
km 2 surface size under water/ice detector. In Section IV, we summarize our results. 



II. INTRINSIC COSMIC TAU NEUTRINO PRODUCTION AND VACUUM 

FLAVOR OSCILLATIONS 

AGNs are the brightest objects in the sky and high-energy photons reaching tens of thou- 
sands of GeV have been observed from them. This is commonly interpreted as an indication 
that some kind of Fermi acceleration is taking place. In conventional models, electrons are 
the particles that get accelerated. It has been argued that if Fermi mechanisms are able to 
accelerate electrons in these objects, protons could also be accelerated by them. In proton 
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acceleration models, the photons arise from neutral pion decays either in jry or pp collisions. 
If this is true, electron and muon neutrinos are also expected to be produced at similar flux 
levels from charged pion decays. Neutrino detection can provide the signature for proton 
acceleration in AGN and it is one of the main goals of neutrino detectors in construction 
and design stages. For an update review of various u e and z/ M flux estimates from AGNs 
in p7 and pp collisions (as well as from some other interesting astrophysical/cosmological 
sites), see jnj. We will here explore the expected tau neutrino fluxes intrinsically produced 
in the collisions and how these fluxes can vary in a possible neutrino oscillation scenario. 

In p7 collisions, high-energy v e and v ^ are mainly produced through the resonant reaction 
p + 7 — > A + — > n-\-7i + . The same collisions will give rise to a greatly suppressed high-energy 
v T (and v T ) flux mainly through the reaction p + 7 — > Dg + A + D°. The production 
cross-section for Dg is essentially up to three orders of magnitude lower than that of A + 
production for the relevant center of mass energy scale. Moreover the branching ratio of 
to decay eventually into v T (p T ) is approximately two orders of magnitude lower than for A + 
to subsequently decay into u e and through n + . These two suppression factors along with 
the relevant kinematic limits give approximately the ratio of intrinsic fluxes of tau neutrinos 
and muon neutrinos as: (u T + v T )/{v^ + u^) < 10~ 5 . 

In pp collisions, the v T flux may be obtained through p + p — > Dg + X. The relatively 
small cross-section for Dg production together with the low branching ratio into v T implies 
that the v T flux in pp collisions is also suppressed up to 4—5 orders of magnitude relative to v e 
and/or fluxes. The situation is quite similar to the prompt atmospheric v T flux calculation; 
the result is basically a rescaling of the prompt flux from the decay of charmed D's and 
results in a negligibly small v T flux for the energies under discussion |T5|,|IE|] . 

In proton acceleration models the intrinsically produced tau neutrino flux is thus expected 
to be very small, typically a factor between 10~ 5 and 10 -6 relative to electron and muon 
neutrino fluxes [lTj . However, recent experimental measurements of atmospheric neutrinos 
suggest that neutrinos could just have vacuum flavor oscillations and the tau neutrino flux 
would be dramatically enhanced. 

It has been pointed out that there are no matter effects for high-energy cosmic tau 
neutrinos originating from cores of AGNs primarily because of relatively small matter density 
in the vicinity of core of the AGN for all relevant 5m 2 [pl|. We will restrict the following 
discussion to vacuum oscillations between two flavors, and u T , for simplicity. 

The flavor precession probability for non vanishing vacuum mixing angle is obtained from 
the effective Hamiltonian matrix in the two flavor basis ip T = (v^, u T ): 



( f sin 26 ^ 
1 f sin2fl 5cos29 , 



(1) 



where 5 = 5m /2E with 5m = m r(y T ) — to and E the neutrino energy, leading to the 
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well known result: 

P{y^v T ) = sm 2 29sm 2 (^l)j. (2) 

If we take the values of sin 2 29 and Sm 2 suggested by recent superkamiokande data (sin 2 29 ~ 
1, 5m 2 ~ 10~ 3 eV 2 ) |TJ| and L ~ 100 Mpc (1 pc ~ 3 • 10 18 cm) as a representative distance 
between the AGN and our galaxy, then the above rapidly oscillating probability averages 
out to ~ 1/2 for all relevant neutrino energies to be considered for detection. Very similar 
fluxes of muon and tau neutrinos would thus be expected. Let us further note that after 
averaging the P given by Eq. (2) is independent of not only E but also Sm 2 and thus leads 
to a constant suppression of high-energy cosmic muon neutrino flux. 

The deficit measured by superkamiokande in atmospheric muon neutrino flux may cur- 
rently be explained either through z/ M — > v T or through — > u a , where u s is a sterile 
neutrinoQ In the first case and for high-energy neutrinos originating at cosmological dis- 
tances, the ratio [y T + z/ T )/(z/ A1 + z/^) is close to 1/2. Therefore, a ratio different from 1/2 
excludes this possibility. 



III. DETECTION OF HIGH-ENERGY COSMIC TAU NEUTRINOS 

High-energy cosmic tau neutrino detection could be achieved by making use of the char- 



acteristic double shower events [|TUJ or by the pileup effect expected as they travel through 
the Earth II] . Such events could be seen in conventional neutrino telescopes and in prin- 
ciple also with other alternative techniques that have been proposed. We will discuss in 
some detail the possibility of detecting double shower events for conventional underground 
telescopes by estimating rates using the the fluxes of Ref. |n| and the oscillation probability 
addressed in the previous section. This is intended to provide a reference calculation. 

The downgoing cosmic tau neutrinos reaching close to the surface of the detector may 
undergo a charged current deep inelastic scattering with nuclei inside / near the detector and 
produce a tau lepton in addition to a hadronic shower. This tau lepton traverses a distance, 
on average proportional to its energy, before it decays back into a tau neutrino and a second 
shower most often induced by decaying hadrons. The second shower is expected to carry 
about twice as much energy as the first and such double shower signals are commonly referred 
to as double bangs. As tau leptons are not expected to have further relevant interactions 
(with high-energy loss) in their decay timescale, the two showers should be separated by a 
clean //-like track [[TO]. 



1 Although, — > v s is now being disfavoured 
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We are going to restrict our estimate to downgoing neutrinos as at these energies tau 
neutrinos that go through the Earth interact. Effectively the process of interaction and 
tau decays can be regarded as an energy degradation to the range 10 4 — 10 5 GeV |TT . 
Unfortunately the two shower signature will be difficult to be resolved below ~ 3 ■ 10 5 GeV 
(see below). 

For downgoing cosmic tau neutrinos, the double bang event rate in water/ice is estimated 
using [^T| 

Rate = A J dEP T (E, ^ nin )^, (3) 

where A is the area of the neutrino telescope and P T gives the probability that a tau neutrino 
of energy E produces two contained and separable showers with the tau lepton energy greater 
than E™ m . It is given by: 



P T (E,Er i ) = Pm -N A dy [D R T ] (4) 

Jo dy 

where Na is the Avogadro's number, p m is the density of the detector medium, da cc /dy is 
the charged current v T N differential cross-section, y is the fraction of neutrino energy that 
is transferred to the hadron in the laboratory frame. D is the detector length scale which 
we fix to be 1 km and the tau lepton range R T which must be contained within D is given 
by: 

R, = (5) 

In Eq. (5), rc is the lifetime and m T c 2 is the mass of the high-energy tau lepton. 

The lower limit of integration in Eq. (3) is E™ lTi . We take E™ m greater than ~ 2 • 10 6 
GeV because at this energy the tau lepton range (separation between the two showers) in 
water is ~ 100 m which allows a clear separation between the two showers. The upper limit 
of integration in Eq. (3) is taken to be E % 2 • 10 7 GeV as for energies above it the tau lepton 
range exceeds the telescope size (see Fig. 1). 

Finally in Eq. (3), dN/ dE is the differential high-energy tau neutrino flux and is obtained 
by multiplying Piv^ — > v T ) given by Eq. (2) with dN /dE for [v^ + u^) taken from fli~3" |. In 
Fig. 2, we depict downgoing differential event rates for double shower events using the parton 
distributions MRS Ri from |22] for km 2 under water neutrino telescopes as an example. We 
have checked that other modern parton distributions give quite similar events rates and are 
therefore not depicted here. We have taken into account the fact that ~ 15% of the times 
the tau decay does not induce any shower. For comparison we also plot the /i-like event 
rate induced by muon neutrinos. Note that these 15% and the tau neutrino interactions in 
which the tau lepton decays outside the detector volume have identical (/z-like) experimental 
signature. 
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The signature of double shower events depends on the detector capabilities for shower 
identification and energy resolution and difficulties can be envisaged. We have used 100 m 
as the minimum distance to resolve two showers, what is quite conservative in view of typical 
spacing between Optical Modules in an under water/ice detector. In Fig. 3, we show the 
dependence of shower size and shower separation on neutrino energy E in ice and in water 
for which we have used the parametrization of ||23| . As shower size is basically proportional 
to energy, the size of the second shower is on average a factor of 2 higher than the first 
one (see also [[TO]] ) . This value results by taking into account the relevant kinematics of the 
allowed decay channels and the corresponding branching ratios and using the average energy 
transfer (y) = 0.25. The y distribution and decay kinematics will lead to a spread in this 
ratio. While y — 0.1 enhances the energy ratio of the second and first showers to a value 
of about 6, for y values higher than 0.4 the ratio of the two shower sizes starts to be lower 
than unity obscuring the tau neutrino signature. 

Another relevant point is the evaluation of the backgrounds, a double shower signature 
not induced by a tau neutrino. As it was discussed in |10j such probability is very small 
and should not affect the detection of the high-energy cosmic tau neutrino. Also one should 
take into account the possibility that the muon component of a single cascade induced from 
a muon or an electron neutrino charged/neutral current interaction can be confused with 
the second shower of the tau lepton decay. However, in this case the size of the second 
shower is smaller than the first one which should be sufficient to distinguish it from a tau 
neutrino event. Thus, the selection criteria of amplitude of second shower greater than the 
first one typically by a factor of ~ 2, depending on y value, essentially makes the observation 
background free. 

The high-energy neutrino telescopes have quite small double shower event rates 
(yr _1 sr _1 ) due to small high-energy intrinsic tau neutrinos flux, thus any observed change 
in this situation may provide indirect evidence of neutrino mass. A corresponding com- 
parable change in this situation is currently not expected from variations of astrophysical 
model inputs. The almost simultaneous measurement of the two showers may provide useful 
information on the incident neutrino energy as well as the y distribution. 

Summarizing, in the context of relevant backgrounds, we envisage essentially the simul- 
taneous presence of two types of events (with different topologies) serving as background for 
the tau neutrino induced contained but separable double showers connected by a /z-like track 
such that the amplitude of the second shower is typically 2 times the first one. The first type 
of background events are due to relatively long (~ 10 km) range muons passing through the 
detector identified as //-like tracks. Their estimated number is given by the upper slanted 
curve in Fig. 2. The second type of background is the single showers due to charged/neutral 
current interactions. These may be estimated as 1/10 of the continuous //-like tracks. Thus, 
the signature of the tau neutrinos as emphasized earlier remained distinct from these two 
type of backgrounds. 
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We emphasize that within the respective energy window, the essential factor in prospec- 
tive detection of contained but separable double shower events connected by a /i-like track 
as a signature of tau neutrinos is the difference in the incident tau neutrino energy depen- 
dences on spread and separation of the two showers. This difference is also clearly crucial 
for separating the tau neutrino events from the (relatively abundent) /i-like events. 



IV. RESULTS 

The intrinsic fluxes of the high-energy cosmic neutrinos originating from proton acceler- 
ation in cores of AGNs are estimated to have typically the following ratios: {y T + z/ T )/(z/^ + 
Vff) < 10~ 5 . Thus, if an enhanced [y T + v T )/{u^ + z/^) ratio (as compared to no precession 
situation) is observed correlated to the direction of source for high-energy cosmic neutrinos, 
then it may be an evidence for vacuum flavor oscillations of neutrinos induced by non zero 
vacuum mixing angle depending on the finer details of the relevant high-energy cosmic neu- 
trino spectra. For vacuum flavor oscillations of high-energy cosmic neutrinos, the relevant 
range of neutrino mixing parameters are: 5m 2 ~ 10~ 3 eV 2 with sin 2 29 ~ 1. 

We have identified the incident tau neutrino energy range and the relevant neutrino mix- 
ing parameters which may give rise to high-energy cosmic tau neutrino induced downward 
contained but separable double shower events. For 2 ■ 10 6 < E/GeV < 2 • 10 7 , a km 2 detector 
may be able to obtain first examples of downgoing high-energy cosmic tau neutrinos through 
contained but separable double shower events or may at least provide some useful relevant 
upper limits. 
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FIG. 1. Comparison of the tau lepton range and the shower length of the first shower (defined 
as twice the depth at maximum) for ice/water. 
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FIG. 2. Expected downgoing /[/-like and r-like event rate produced by AGN neutrinos from ref. 
[13] because of vacuum flavor oscillations as a function of the minimum energy of the corresponding 
lepton. 
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FIG. 3. Typical longitudinal development of a double shower produced by the deep inelastic 
charged current interaction of the high-energy cosmic tau neutrino. 
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